Generalized transformation optics from triple spacetime metamaterials 
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In this paper, various extensions of the design strategy for transformation media are proposed. We 
show that it is possible to assign different transformed spaces to the field strength tensor (electric field 
and magnetic induction) and to the excitation tensor (displacement field and magnetic field), resp. 
In this way, several limitations of standard transformation media can be overcome. In particular, 
it is possible to provide a geometric interpretation of non-reciprocal as well as indefinite materials. 
We show that these transformations can be complemented by a continuous version of electric- 
magnetic duality and comment on the relation to the complementary approach of field-transforming 
metamaterials. 



I. INTRODUCTION 

In the field of metamaterials, artificial electromagnetic 
materials, the use of spacetime transformations as a de- 
sign tool for new materials has been proved very success- 
ful recently [ll, [2, Q ■ As basic idea of this concept a meta- 
material mimics a transformed, but empty space. The 
light rays follow the trajectories according to Fermat's 
principle in this transformed (electromagnetic) space in- 
stead of laboratory space. This allows one to design in an 
efficient way materials with various characteristics such 
as invisibility cloaks [l|, 0, 0] > perfect lenses \3\ , magnifi- 
cation devices Q , an optical analogue of the Aharonov- 
Bohm effect or even artificial black holes Q. Still the 
media relations accessible in this way are rather limited, 
in particular non-reciprocal or indefinite media (materi- 
als exhibiting strong anisotropy) are not covered. But 
these types of materials also have been linked to some of 
the mentioned concepts, in particular perfect lenses [a, |7[ 
and hyperlenses Q. This raises the question whether 
there exists an extension of the concept of transforma- 
tion media such as to cover those materials as well and 
to provide a geometric interpretation thereof. 

In this paper we propose an extension of this type. 
As in Refs. [ESS] our concept is based on diffeomor- 
phisms locally represented as coordinate transformations. 
Therefore many of our result allow a geometric interpre- 
tation similar to the one of Refs. [3|, |9| as opposed to 
another recently suggested route to overcome the restric- 
tions of diffeomorphism transforming media [lO, [ll| . The 
starting point of our considerations are Maxwell's equa- 
tions in possibly curved, but vacuous space [24| : 

V.,B' = , VoB' + e'l'^djEk = , (1) 

V{D' = p , e'^'^djUk ~ VoP^ = f . (2) 

Here, Vi is the covariant derivative in three dimensions 

V.A^ = (9. + T\^)A= = -^a,(V7^^) , (3) 



with the space metric 7^^ and its determinant 7. 

For many manipulations it will be advantageous to use 
relativistically covariant quantities. Therefore, Eqs. ([1]) 
and ([2) are rewritten in terms of the field strength tensor 
F^^, the excitation tensor iJ^"^ and a four current J^ (cf. 
Appendix ^ : 



e^^'^P'^d^F.^^O , D,H^"' 



-JP 



D^J^ = 



(4) 



The four dimensional covariant derivative D^ is defined 
analogously to ([31), whereby the space metric is replaced 
by the spacetime metric g^i, and its volume element \J—g. 
We wish to analyze these equations of motion from 
the point of view of transformation media. All transfor- 
mation materials have in common that they follow as a 
transformation from a (not necessarily source-free) vac- 
uum solution of the equations of motion, which maps 
this solution onto a solution of the equations of motion 
of the transformation material 25]. The crucial ingredi- 
ent in the definition of transformation media then is the 
class of transformations to be considered. As space of all 
transformations we restrict ourselves to all linear trans- 
formations in four-dimensional spacetime. Consequently, 
all media exhibit linear constitutive relations, which may 
be written within the covariant formulation as Il2l 






In vacuum one obtains l2f 



f 



= 7T (.9'"'5''" - g'^^g""') 



(5) 



(6) 
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such that the standard result E — f) and B = H emerges. 
These transformations and the ensuing media proper- 
ties ([5]) have the advantage of being relativistically in- 
variant and thus very easy to handle. However, they do 
not include any frequency dependence and remain strictly 
real, which perhaps is the most severe restriction that fol- 
lows from the coordinate transformation approach. As 
long as the linear transformations are seen as transfor- 
mations of spacetime (rather than of the fields) this re- 
striction is not surprising, though. Indeed, from energy 
conservation it follows that it is impossible to model a 



process of absorption by the medium as a local transfor- 
mation of spacetime (notice that the spacetime itself is 
not dynamical and thus cannot contribute to the energy) . 



\e, b] {x) diffco- 



\e,b]{x) 



\E,B\(i = x(x)) 



KD,H (x) morphism I),W (x) pretation \t> ,n\ (i = x{x)) 



II. DIFFEOMORPHISM TRANSFORMING 
METAMATERIALS 

Obviously, the concept of transformation materials as 
sketched above is related to symmetry transformations, 
as those are by definition linear transformations that map 
a solution of the equations of motion onto another one. 
Therefore it is worth working out this relation in some 
more detail. 

A symmetry is a transformation which leaves the 
source-free [23| action of the theory, here 



S 



rji^u 



(7) 



invariant, whereby surface terms are dropped. It 
straightforwardly follows that a symmetry transforma- 
tion applied on a solution of the equations of motion still 
solves the latter. In the above action a general, not nec- 
essarily flat, spacetime is considered. The symmetries of 
this action are well known: these are the [/(I) gauge 
symmetry of electromagnetism and the symmetries of 
spacetime (diffeomorphisms) . The gauge symmetry can- 
not help in designing materials as the media relations 
are formulated exclusively in terms of gauge invariant 
quantities. However, diffeomorphisms change the media 
relations, as is pointed out e.g. in Ref. flj] and as it has 
been applied to metamaterials in Ref. [3|. Thus one way 
to define transformation media is: 

Definition 1. A transformation material follows from 
a symmetry transformation applied to a vacuum solution 
of Maxwell's equations. This vacuum solution need not 
he source free. 

The space of all possible transformation materials of 
this kind has been derived in Ref. [3|; here we briefly 
want to summarize the result of that paper. The starting 
point is the observation that a curved space in Maxwell's 
equations looks like a medium. Indeed, in empty but 
possibly curved space the constitutive relation among the 
electromagnetic fields is found by exploiting 

Fai = {googij - goj9ia)H°^ + gokguH^'^ , (8) 
ff'-'^25V'i^ofe+gV'i^H, (9) 

which in terms of the space vectors reads 



V^9ao ffoo 



B' ^ 



V~9oo goo 



(10) 
(11) 



Thus empty space can appear like a medium with per- 
meability and permittivity e'-* = /i*^ = .9*"'/V~.9oo a-nd 
with bi-anisotropic couplings ^^-^ = — k*-' = e^^^ goi/goo- 



FIG. 1: Diffeomorphism transforming metamaterials accord- 
ing to Ref. |3|. 



Now, as the basic idea of Ref. Q , if empty space can 
appear like a medium, a medium should also be able to 
appear as empty space. One starts with electrodynam- 
ics in vacuo, we call these fields F^i, and H^^'^' with flat 
metric g^i,. Now we apply a diffeomorphism, locally rep- 
resented as a coordinate transformation x^ — + x^^ {x) . As 
the equations of motion by definition are invariant un- 
der diffeomorphisms, all relations remain the same with 
the fields F, H and the metric g^^ replaced by the new 
barred quantities. As a last step one re-interprets in the 
dynamical equations ([T|) and ^ the coordinates x^^ as 
the original ones x^ , while keeping g^n in the constitu- 
tive relation. To make this possible some fields must be 
rescaled in order to transform barred covariant deriva- 
tives (containing g) into unbarred ones (containing g.) 
The situation of diffeomorphism transforming metama- 
terials is illustrated in Figure [H which also summarizes 
our notation. As a more technical remark it should be 
noted that this manipulation is possible as we consider 
just Maxwell's theory on a curved background rather 
than Einstein-Maxwell theory (general relativity coupled 
to electrodynamics.) In the former case the metric is an 
external parameter and thus this manipulation is possi- 
ble as long as none of the involved quantities depends 
explicitly on the metric. 

To keep the whole discussion fully covariant the fields 
are transformed at the level of the field strength and ex- 
citation tensor (rather than at the level of space vec- 
tors as was done in Ref. 3] ) . To transform the covariant 
derivatives Z)^, into the original D^ we have to apply the 
rescalings 



H^"" 



-.H"" 



Jf" 



-.J^ 



(12) 



In addition the transformation x^ — > a;'' (x) may not pre- 
serve the orientation of the manifold, which technically 
means that the Levi-Civita tensor changes sign [S^]. This 
is corrected by introducing the sign ambiguity 



F„, 



±F„ 



(13) 



with the plus sign for orientation preserving, the mi- 
nus for non-preserving transformations. These new fields 
again live in the original space with metric 5^^, but now 
the space is filled with a medium with 



±-- 



(r^r 



g'-^g"") , 



(14) 



or, in terms of space vectors, 



V' 



V-500 ^/7 



E, - ^e^^'Ui , (15) 



500 



V-900 V7 500 



(16) 



with s = ±1 being the sign in (fT5|) and ([H)) . As can 
be seen, the media properties are restricted to reciprocal 
materials (e = e^ , fj, = fj?" , k = x^)^ which, in addi- 
tion, obey e = fj,. This result has been obtained in Ref. 
Q in a slightly different way and encompasses the trans- 
formations in Refs. [l|, [ij]. We do not want to go into 
further details of this approach but refer to the review 
Q , where its geometric optics interpretation is discussed 
in detail. Indeed, light travels in transformation media 
of this type along null geodesies of the electromagnetic 
space S^, which allows (with some restrictions to be dis- 
cussed in Section|V| a simple and intuitive interpretation 
of the transformation. 



III. TRIPLE SPACETIME METAMATERIALS 

Despite the variety of applications of diffeomorphism 
transforming metamaterials some results suggest a search 
for extensions. Indeed, there exist e.g. designs of super- 
and hyperlenses that make use of indefinite materials 
(strong anisotropy) [a, 0, Q- Though both concepts 
should be perfectly understandable in terms of trans- 
formation media, the specific material relations used in 
these works do not fall under the class of diffeomorphism 
transforming metamaterials. 

To understand a possible route to generalize the con- 
cept of diffeomorphism transforming media we have to 
consider again their basis, namely symmetry transfor- 
mation. The concept of symmetries is used to identify 
different solutions of the equations of motion that effec- 
tively describe the same physics. By means of the re- 
interpretation in the last step of Figure [1] such symme- 
try transformations can be used as a simple tool to derive 
within a restricted class of constitutive relations new me- 
dia properties in a geometrically intuitive and completely 
algebraic way. 

Nonetheless, within the concept of metamaterials it is 
not important that the transformed solution in principle 
describes the same physics as the original one. Still, one 
may want to keep the possibility of mapping source free 
solutions onto other source free solutions in a straight- 
forward way, as only in this way do we have an effective 
control over passive media and do not risk introducing ex- 
otic sources such as magnetic monopoles. Furthermore a 
geometric interpretation of the transformations is kept, 
which is advantageous in many applications. To weaken 
the conditions on transformation materials while keep- 
ing the advantages of symmetry transformations we thus 
propose the following definition: 



Definition 2. Consider the set of all transformations 
T which map a source free solution of the equations of 
motion ([1]) and ([2]) onto another source free solution. A 
transformation material is a material obtained by apply- 
ing a transformation T onto a (not necessarily source 
free) vacuum solution. 

There are two types of extensions contained in this 
definition compared to the previous section: 

1. There exist transformations that leave the equa- 
tions of motion invariant, but change the action by 
a constant and thus are not symmetry transforma- 
tions. A transformation of this type is the so-called 
electric-magnetic duality. Its effect will briefly be 
discussed in Section ITll Al 

2. We do allow for transformations which leave all 
Maxwell's equations ([T]) and ^ invariant, but 
change the media relations fS]). This indeed gen- 
eralizes the concept in an important way. 

To see the origin of the second extension it is impor- 
tant to realize that the equations of motion of electro- 
dynamics separate into two different sets (Eqs. ((T|) and 
([2]), resp.) with mutually exclusive field content. This 
characteristic is not just an effect of our notation, but 
as has been shown e.g. in Refs. [1^1 13 j the equations of 
motion of electrodynamics can be derived from first prin- 
ciples without using explicitly the constitutive relation 
H — H{F). As the two sets of equations are separately 
invariant under diffeomorphisms it should be possible to 
assign different transformed spaces to H — (T>, Ti.) and 
F = {E, B). In other words, it must be possible to distort 
the spaces (or the coordinates) of the field strength tensor 
and the excitation tensor separately, whereby the result- 
ing transformation material per constructionem satisfies 
all conditions of the Definition [51 The ensuing constitu- 
tive relation as well as the solutions of the equations of 
motion still follow (almost) as simple as in the case of 
Ref. [|. 

To prove the potential of this method we have to ex- 
tend the notation compared to the previous section: as 
before laboratory space has metric g^j/ , its fields in vacuo 
are H — {V, H) and F — {E, B); the fields of the trans- 
formation material (living in the space with metric g^^) 
are again labeled with a tilde. The transformed space 
of the field strength tensor has metric g^j/ and fields 

F = {E,B), the one of the excitation tensor ^^jy and 

H = {T>,'H). This new transformation is illustrated in 
Figured Applying the two transformations 



x^' ^x^'{x) , f" ^x^'ix) 



(17) 



to the constitutive relation ([5]) with x being the vacuum 
relation ^ yields 



2 dx^ dx-" 



(18) 



911" 

[e,b]{x) 

\v,n] (x) 




\e,b\{x) 



9;.^ 



[v, h] (x) <f 




that 



E,b] {S: = x{x)) 
V,n\ (i = f(a;)) 



FIG. 2: Illustration and notation of the generalized "triple 
spacetime metamaterials" . Notice that the diffeomorphism I 
acts only on the fields E and B, while diff'eomorphism II acts 
on D and TC. 






<t^ 



s^- 



v~5 



^e,u{c^</=-!p^'), 



2 x/7 
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( ^fik „il „Ol„ik\ 

<jki[g g -g g ) 



(24) 
(25) 
(26) 
(27) 



which are the defining tensors of the Boys-Post relation. 
After some algebra the Tellegen relation 



Introducing the notation 



^ dxPdx'"'^ ^ dxP dxP^ 



a = n ~ a^nd (l-^J 



the relation may be written as 



H''" = 2 (^^"5'* " ^''*5'') ^"^ 



(20) 



It should be noted that g^^"^ in Eq. (I19p is no longer a met- 
ric, in particular it need not be symmetric in its indices 
and it need not have signature (3, 1). 

To derive the new constitutive relations in the original 
(laboratory) space we proceed analogously to the previ- 
ous section. All fields have to be rescaled in order to obey 
the equations of motion in the original space with metric 
5^y, which implies 



\/75oo 75oo 



3 ' 

(28) 



B^ = -S^^^cf^j + s-|^^^^-^/-efei™g'o/^^^£:, 



V^m" 



19m 



(29) 



is found, which in the limit of g^jy — g^^, is equivalent to 
Eqs. (jlSp and (|16p . An important comment is in order: 
due to the different transformations applied to RP" and 
Ffj^jy, resp., the constitutive relation ((22)l . or ((28)) and (|29p . 
relates fields from different spacetime points in the orig- 
inal space, e.g. Ei {x = x{x)) refers the field Ei{x) at a 
different point x^ in the original space than I?* {x = x{x)) 
does. 

Let us comment on the more technical parts of this 
result. In Section[TT]we saw that transformation materials 
derived from symmetry transformations are restricted to 
reciprocal materials with e = fi. These restrictions can 
be overcome partially with the above result: 



F^, = ±F^, , iJ'"- = X^H^''' , > 



Thus the constitutive relation becomes 



^J^ 



(21) 



HP-^X^^-P-Fp„^±-- 



-g 



ig'''g"'~g'"'g''lFp., 



where the sign refers to the possible change of orienta- 
tion in the transformation x^ ^^ x^. For the equivalent 
relation in terms of space vectors the notation 



(23) 



is used, where s and s are the respective signs due to 
the change of orientation in the transformations to lab- 
oratory space. Now it easily follows from (|A.17p - (|A.20p 



• As g*-' = [g^'^y it follows that permittivity and 
permeability are related as 



5M'^ = s\ 



-ge 



ji 



(30) 



It should not come as a surprise that permittiv- 
ity and permeability cannot be independent, as by 
virtue of the definition of the relativistically covari- 
ant tensors i^^j/ and if^i/ such transformations can- 
not act independently on E and B or V and Ti, 
resp. A possible route to relax this restriction is 
discussed in Section IVll 



• Permittivity and permeability need no longer be 
symmetric. Therefore it is possible to describe 
non-reciprocal materials, or, in the language of Eq. 
(JA.23P , the skewon part need not vanish. This hap- 
pens if the mapping between the two electromag- 
netic spaces, dxP'/dx'^, is not symmetric in fi and 
v, e.g. for a material with mapping x = x — z, 
x = x + z. 



The generalized transformations yield many more 
possibilities considering the signs of the eigenval- 
ues of permittivity and permeability. Within the 
method of Ref. [3], /i and e are essentially deter- 
mined by the spatial metric of the electromagnetic 
space (cf. Eqs. (fT5|l and (fT6|) and recall the relation 
gij — ^*i.) However, a spatial metric by definition 
must have three positive eigenvalues, a characteris- 
tic that cannot be changed by any diffeomorphism. 
Thus it follows that in any medium of this type the 
eigenvalues of e and /j, are all of the same sign. 

— Within the generalized setup of "triple space- 
time metamaterials" , however, the signs of the 
eigenvalues in e can be chosen freely, as the 
metric is multiplied by a transformation ma- 
trix, 






(31) 



and no restrictions on the signs of the eigen- 
values of the transformation matrix exist. In 
this way indefinite materials [g, 0| can be de- 
signed as a result of different space inversions 
in the two different mappings. As an ex- 
ample the mapping z = — z, z — z (with 
all other directions mapped trivially) yields 
e^^" = diag(-l, -1, 1), ^i'^ = diag(l, 1, -1). 

Furthermore the relative sign between the 
eigenvalues of e and those of /i can be chosen 
as a consequence of the factor s in Eq. ([30)1 . 
This change in the relative sign may be inter- 
preted as a partial reversal of time as can be 
seen in the following list (space maps trivially 
here and all media are assumed to be homo- 
geneous) : 
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-1 



We note that all eight classes of materials discussed 
in Ref. [g| allow a geometric interpretation within 
the setup of "triple spacetime metamaterials." 

• More complicated than permittivity and perme- 
ability are the bi-anisotropic couplings. With the 
standard assumption of g'^^ = in laboratory space 
it follows from 



-ss- 
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g g ekimg 
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(32) 
(33) 



similarly to Eq. ([151) that all electric-magnetic cou- 
plings vanish if the transformation does not mix 
space and time. In this case the crucial compo- 
nents g°' and 5°' may be written as 



dx'^ 00 "5^ 



9 



dx° ^j dx'- 
dx° dx^ dx^ dx^ 
9x° 00 ^^ 



dx° ^j dx'- 



(34) 
(35) 



Most importantly it is found from these expres- 
sions that one of the two bi-anisotropic couplings 
may vanish while the other one is non- vanishing, 
which is impossible within the context of diffeomor- 
phism transforming media. Moreover, in the latter 
case the bi-anisotropic couplings must be symmet- 
ric matrices, which need no longer be the case in 
the present context. 



Finally, the result (pS)) and ^^ reduces to the rela- 
tions (J15p . (|16p if gpp is a symmetric matrix of sig- 
nature (3, 1) and, in addition, \/—g — \J —g. This 
does not necessarily imply x^ = x^ but rather that 
there exists yet a different space which describes 
the same media properties in terms of the transfor- 
mations of Section jTTI. 



A. Electric-magnetic duality and rotation 



Finally, we should ask whether Eqs. ([28)) . ([29|1 indeed 
describe the most general media fulfilling Definition |21 
Taken separately, the two sets of equations in (|T]) and ([2]) 
do not exhibit more symmetries than diffeomorphisms. 
However, there exists the possibility of transformations 
that mix F^^^ and i?'"'. Indeed a transformation of this 
type is known as electric-magnetic duality, which has im- 
portant implications in modern theoretical high-energy 
physics [I3|. It represents the fact that under the ex- 
change 



I iiv ■*"* *Iii.ii' 



or in terms of space vectors 






'Hi , 



n, 
v- 



'E,, 



B' 



(36) 

(37) 
(38) 



the source-free equations of motion do not change (the 
action changes by an overall sign.) Of course, this duality 
transformation is problematic when applied to a solution 
with sources, as it transforms electric charges and cur- 
rents into magnetic charges and currents and vice versa. 
In the remainder of this section we thus restrict to source- 
free solutions or should allow the possibility of artificial 
magnetic monopoles. Then it can be checked straight- 
forwardly that electric-magnetic duality applied to the 
result ([^01) . or ([^5)1 and ([^^ . does not yield media rela- 
tions not yet covered by diffeomorphisms alone. 
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FIG. 3: Illustration and notation of the generalized "triple 
spacetime metamaterials" complemented by electric-magnetic 
rotation. The electric-magnetic rotation must act after the 
transformation of spacetime as these two steps do not com- 
mute. 



However, as far as the equations of motion ([T]) and 
are concerned, electric-magnetic duality can be promoted 
to a continuous U{1) symmetry with transformation '28|| 



IV. PERFECT LENS FROM INDEFINITE 
MATERIAL: AN EXAMPLE 

To provide a better understanding of the formalism 
developed in the previous section a concrete example is 
demonstrated. To keep things simple we show how a pro- 
posal taken from the literature can be given a geometric 
interpretation. 

In Ref. [H it has been pointed out that two slabs of in- 
definite material (media with strong anisotropy) can form 
a perfect lens. Since, in constrast to standard diffeomor- 
phism transforming media, strong anisotropy is available 
in triple spacetime metamaterials the question appears 
whether a geometric interpretation of the lens proposed 
in Ref. [(I can be given (cf. Ref. [I^ for a related discus- 
sion.) We consider the lens to be an infinite slab in the 
x-y plane with a certain thickness in the z direction. In 
its simplest form the lens consists of two slabs of equal 
thickness d, where the media properties of the first slab 



B"" = cos aB' - sin aV , V = cos aV + sin aS* , 

(39) 

Ei = cos aEi + sin aJii , Tii = cos oHi — sin aEi . 

(40) 

These transformations comply with Definition [2] and thus 
their action onto a medium with general constitutive re- 
lation (jA.13|) should be studied. The result 

V^ = (cos^ ae + sin^ afi + sin a cos a{K + £,)) Ej 

-f- (cos aK — sin a^ -f- sinacosa(/i — e)) Hj , 

(41) 

i?* — (cos^ a/1 + sin^ ae — sin a cos a{K + C)) "^j 

+ (cos^ a^ — sin^ an + sin a cos a(/x — e)) Ej 



(42) 



shows that the transformation acts trivially if e — fi and 
^ — — K, in particular in vacuo and consequently for all 
diffeomorphism transforming media (jlSp . However, they 
yield new media relations when acting on a solution of 
the type (pS)) and (P^ . Therefore these new relations 
are part of the materials covered by Definition [2l They 
are derived here for completeness, though their geomet- 
ric interpretation is not immediate. The coordinate lines 
x'^(x) and x'^{x) could be understood as the electromag- 
netic spaces of the linear combinations {E, B) and (2?, 7Y) 
as given in ([M)l and (PD|) . resp. Still, one should be careful 
with this interpretation: as the transformation of space- 
time does not commute with the electric-magnetic rota- 
tion one cannot modify the situation in Figure[3]in such a 
way that the two electromagnetic spaces, x^ and S^, are 
identified with certain linear combinations of {E, B) and 
(I?,7i), resp.; rather the electric-magnetic rotation acts 
upon the fields after the transformation of spacetime. 



e'^=A**-'=diag(l,l,-l) 
while in the second slab 



e'-'=M^^=diag(-l,-l,l) 



(43) 



(44) 



To provide a geometric interpretation we start with the 
observation that a standard perfect lens with e = /i = — 1 
may be produced by two different transformations, either 
a space inversion £ = — z, or a time reversal t = —t. 
From Eqs. P^ and (fTB|) it follows straightforwardly that 
these two transformations yield the same media proper- 
ties. Within triple spacetime metamaterials we now may 
ask the question of what happens if space inversion is 
applied to one set of the fields, while time reversal is ap- 
plied to the other set. For concreteness, space inversion 
is applied to the fields E and B and thus 



z = —z -|- Z\ 



(45) 



where Z\ is an unimportant constant necessary to meet 
the boundary conditions. All other fields x^ are mapped 
trivially. The second set of fields, T) and 7Y, transform 
according to 



i = -i + Ti 



(46) 



with all other fields transformed trivially. Consider now 
these two transformations in Eqs. (P5| and (P5)) . From 
(irai) one finds 



7'^=.9'^=diag(l,l,-l) 



(47) 



Furthermore, s = s = — 1 as both tranformations are 
orientation changing. Furthermore, f/gg = 1 (remember 
our convention ^oo = ~l)i such that indeed the media 
properties (|43p are found in this slab. 

This single slab of indefinite material does not establish 
a perfect lens, as can be seen easily when studying how a 




FIG. 4: Mapping of the world-line s{t) = (t, z{t)) in the origi- 
nal space onto the deformed spaces by means of the two differ- 
ent transformations. Blue lines indicate the trajectory outside 
of the lens (trivial mapping of all x'^), the red line represents 
the first slab, the green line the second slab. The parametriza- 
tion of the world-line in the original space is assumed to obey 
z{t) =t. 



world-line s{t) = (t(T),a;(r), y(T), z(t)) is mapped onto 
the two deformed spaces. For simplicity time may be 
interpreted with the parametrization variable ^(t) = r 
and furthermore we can assume without loss of generality 
z{t) = t = t [2§|. The situation is illustrated in Figure ID 
As can be seen the mappings do not agree after the first 
slab, both trajectories are at the same point z = 0, but 
they differ in time. This must be corrected in the second 
slab. In our example we have chosen a completely trivial 
mapping for V and Ti, so these fields propagate in the 
second slab as in free space. E and B, however, are 
transformed as 



t 



-t + n 



-z + Zo 



(48) 



which actually reverts the transformation (|45|) and at the 
same time applies (|46|). Not surprisingly, the two trajec- 
tories now meet at the same point again and the perfect 
lens is established. Again it is immediate that this trans- 
formation establishes the media relations (|^^ . There- 
fore, triple spacetime metamaterials indeed can provide 
a geometric interpretation of the lens of Ref. [6| . It should 
be noted, that this specific lens has focal length zero, it 
shrinks the effective width of the device from 2D to zero, 
but not to a negative value as is necessary for a real lens. 



ENERGY, MOMENTUM AND WAVE 
VECTOR 



So far we studied solutions of Maxwell's equations 
which — up to rcscalings — are equivalent to certain vac- 
uum solutions. Still we did not ask up to what point these 
transformation materials really are "media that look like 
empty space." To do so it is not sufficient to consider the 
transformation of the fields and sources, but equally well 
we should look at the conservation laws, summarized in 
the conservation of the stress-energy-momentum tensor 



(SEM tensor.) While in the generic situation of electro- 
dynamics in media, the definition of the "SEM tensor of 
electrodynamics" is not unique [1^, [23| , we do not have to 
deal with these subtleties in the present situation as our 
(idealized) media are lossless and dispersion free and thus 
allow for a definition of a complete action (cf. Eq. ([7])) 
without any reference to "matter." Therefrom we imme- 
diately derive the covariant SEM tensor 



rp^V 



1 

47r 



Fp„g''^^^H''^P + -a^^'^ F.^^Hf' 



pa^ 



D^T'"' = 



(49) 
(50) 



The advantage of this tensor over the canonical SEM ten- 
sor is the simple behavior under diffeomorphisms: being 
a real tensor field, T^i, transforms exactly in the same 
way as the metric. 

Let us now look at the materials as described in Sec- 
tion|TTl Thanks to its transformation properties the SEM 
tensor in the electromagnetic space follows immediately 
as 



An \ 4 



p(jj 



(51) 



But how about T^^7 Of course one could define an "in- 
duced SEM tensor" from the electromagnetic space as 
(cf. Eqs. lEl) and ^) 



fr = T- 



-gAn 



Fp^g^'^H'^^P + -gf^'^Fp^HP" 



(52) 

but obviously this tensor is not conserved in laboratory 
space, DpTj'^ ^ 0, since it depends explicitly on the met- 
ric gfiu- In other words, the crucial trick to re- interpret in 
the dynamical equations the coordinates in electromag- 
netic space, a;^, as those in laboratory space, x^ , works 
in the equations of motion ^ and ([2]) , but does not work 
for the SEM tensor and its conservation. 

Of course, the correct SEM tensor in laboratory space 
immediately follows from (|49p as 



f" = -— (Fp^g''^t'H''^P + -g^^Fp^HP" 
An \ 4 



(53) 



Clearly, requiring equivalence of the two tensors would 
not even allow for conformal transformations. But even 
when looking at integrated quantities (total energy and 
momentum flux in the material), 



P'^ = fd^x^T°'' 



(54) 



the induced tensor does not yield the correct quantity in 
laboratory space. Of course, the situation is even more 
complicated for triple spacetime metamaterials: since the 
transformation of the explicit metrics appearing in Eq. 
(H5)l is not defined, an "induced SEM tensor" cannot even 
be defined. 



Instead of the correct, directly evaluated SEM tensor 
([55)1 a slightly different tensor is considered in the follow- 
ing. To see its advantage we make the standard assump- 
tion that our laboratory space metric has goo = ~1 {x'^ 
is our laboratory time) and goi = (the measure of dis- 
tances is time independent.) Then it is straightforward 
that the quantity 



M'"' 



-g^PFp^H'^'' 



(55) 



contains the Poynting vector and the direction of the 
wave vector, 






(56) 
(57) 



We recall that the original fields obey the constitutive 
relations of vacuous space and thus trivially Af"' = Af". 
From the transformation rules (pij) the transformed ten- 
sor is found as 



M"" 



'.g^'Fp^H" 



-9 



(58) 



^^ dxP dx'' dx" 



H^ 



dxP ■ 



The transformation law of M^'^ encodes in a geometric 
language how energy flux and phase velocity behave in a 
medium. For simplicity let us now concentrate on media 
without bi-anisotropic couplings, in other words we allow 
for general spatial transformations as well as stretchings 
and reversal of time, but keep g^i = goi = 0. Then we 
find for the transformed space vectors (cf. Eqs. (|A.9p - 
(IXT2]) ): 



S" 






-Ss^^^'j^e-^'^E^j^n.^ (59) 
500 ox'^ ox^ ox'^ 






dx" 
I 

dx^ 



^/7 dx^ 



- — P^^B" (60) 



The transformation of the Poynting vector may be ab- 
breviated as 



and it is then easily seen that rii transforms as 
__==\/77 dxP dxP 



where C/,,,i, is the inverse of T*^*^ in the sense of 



Hjk 



U^.kT'^'' 



(61) 



(62) 



(63) 



While these formulae might look cumbersome, their geo- 
metric interpretation actually is quite straightforward. 
In the case of diffeomorphism transforming materials, 
gfj.i' = gtiu, Eq. (|55|) states that S** behaves under purely 



spatial transformations as a covector [l|, while rii from 
Eq. ([60|) behaves as a vector: 



S'^ 



\f~g dx^ dx^ 
\/-9qoV^ dx° dx^ 
, — dx^ dx^ 



S^ 



'^ ^dx^dx^'"' 



(64) 
(65) 



Of course, the relative orientation of 5* and n' is pre- 
served under the diffeomorphisms, but this is no longer 
true for S"* and rV' , since indices are raised/lowered by the 
space metric 7^ in laboratory space as opposed to 7^ in 
electromagnetic space. 

For triple spacetime metamaterials no linear transfor- 
mation 5' — T'^jS^ exists. This makes the interpretation 
a little bit more complicated, but at the same time is the 
source of the numerous additional possibilities within this 
generalized setup. In general, the value of the element 
T'^*^ defines the component of the Poynting vector in di- 
rection x' as generated by electric and magnetic fields 
that point in the original space in the directions x^ and 
x'^ , resp. In this way it is easy to engineer the direction 
of the Poynting vector in the medium for a given polar- 
ization of the incoming wave in vacuum. Similar conclu- 
sions apply for the transformation matrix Uijk , with the 
notable restriction that n* can be parallel or anti-parallel 
to k\ Whether ('D^S-', A:') form a right- or left-handed 
triple can be deduced from 



eyfcP^B^ = ^E.e^'^Ek = 



\/75o5 



'^E.y'E, 



(66) 



A. Wave vector and dispersion relations 

While the above relations correctly reproduce the di- 
rection of the Poynting and the wave vector, they cannot 
distinguish between propagating and evanescent modes. 
Consider as an example the following transformation: 



iXj <*-' • iXj iX' 



=0 _ _ =i _ i 



(67) 



From Eqs. (|T5|) and (|16p it is found that this is a homo- 
geneous material with e — ~1 and /i = 1. As fields in 
vacuo, E — V and B — Ti, we consider a monochromatic 
wave 

E ^ ee*(^^--*) + c.c. , fc • e = , (68) 

B = 6e^('=^-'^*) + c.c. , b=-kxe . (69) 

UJ 

After the transformation the fields E, B and V, Ti refer 
to the original fields at different time instances: 

Eixf" =x^'{x))^E{x,t) , (70) 

i3(i'^ =i^(a;)) ^B{x,t) , (71) 

Vii^" =x''{x))^-V{x,-t) , (72) 

H {i" = xf'ix)) = -nix, -t) . (73) 



Of course, Maxwell's equations are satisfied by the new 
fields by construction. Still, the partial exchange of pos- 
itive and negative angular frequencies has important im- 
plications in the dispersion relation as any propagating 
wave in vacuo becomes evanescent in the medium and 
vice versa. 

Though this behavior may not appear immediate when 
transforming the monochromatic wave (|68p and (|69p with 
([70)1 - ([75]) ■ it can be made explicit from geometric quan- 
tities as well. Indeed, from the relativistic wave equation 

M 



D^x'^'P'^DpA, ^-r 



(74) 



it follows straightforwardly that "triple-space metamate- 
rials" in the absence of charges and currents and in the 
limit of approximate homogeneity obey the dispersion re- 
lation 



7 KnKu — U , 



{LO,k) 



(75) 



In our example the partial reversal of time yields g'"^ = 1 
and thus lo'^ + lP — 0. 



VI. NON-INVARIANT TRANSFORMATIONS 

Within the approaches to transformation media dis- 
cussed so far invariant transformations of the equations 
of motions were used exclusively. This means that the 
transformations "do not introduce charges or currents", 
in other words the transformation medium based on a 
source free vacuum solution will be source free as well. 
What happens if this restriction is abandoned? Still in- 
sisting on a constitutive relation of the form ([5]) this sug- 
gests the following definition: 

Definition 3. A transformation medium is defined by 
an arbitrary linear transformation applied to a (not nec- 
essarily source free) vacuum solution of Maxwell 's equa- 
tions. The linear transformation constitutes the media 
properties as well as charges and currents of the trans- 
formation medium. 

Though not in its most general form, this approach was 
proposed in [IQ. .11.]. Starting from the vacuum relation 
^ the most general linear relation can be achieved by 
the field transformations [301 

H^'' = n^%^HP- , F^, = ^^^.P^Fp^ , (76) 

with the transformed x, 



X'^ 



= (rix*)^'-'' 



(77) 



The original fields F^^^ and H^'^ by assumption are solu- 
tions to the equations of motion. If we allow besides the 
standard electric four-current J^' also a magnetic four- 
current J^j, any transformation of the type (|76p can be 



mapped on a solution of the new equations (^^^'"^ is the 
^A."^ = {5P,5Z - 515P)I2) 



inverse matrix ^ ^^^'^ 



D^H^'' ^.r ^ D^{n^%,HP-) 



(78) 



eP-'P'd^Fp, = J^ = eP-'P-d.i^^p^^'Fr^) , (79) 

provided appropriate currents are introduced. These 
transformations in general are not symmetry transforma- 
tions and accordingly a source free solution is no longer 
mapped on another source free solution. The transfor- 
mations (|76p and the ensuing equations of motion (|78p 
and (|79|) are the relativistic form of the transformations 
proposed in Refs. [10, UM ■ 

We recover the invariant transformations of Section 
pup by introducing the restrictions 



i i per — ^ p^ 



iT; per _ T^ PT ° 



(80) 



Let us first count the degrees of freedom in the transfor- 
mations, x^'''"^ is a rank four tensor, anti-symmetric in 
(/i. u) and (p, a) and thus has 36 independent components 
(20 components of the principal part, 15 of the skewon 
part and one axion coupling). The same applies to the 
transformation matrices Vt and ^. Restriction to diffeo- 
morphisms according to Eq. (|80p reduces this number to 
6 parameters for each, D, and ^; the electric- magnetic 
rotation of Section [ill Al adds another parameter in form 
of a rotation angle. Here, another important difference 
between a transformation material according to Defini- 
tion [21 and [3| emerges. In both cases the transforma- 
tion yielding certain media properties is not unique. In 
the former case, however, different transformations are 
physically equivalent as they are connected by symme- 
try transformations (isometrics of the laboratory metric 
gp,^). In the more general case of Eq. ([7S)) the different 
transformations need not be physically equivalent. As 
is immediate from Eq. (|76p a certain medium exhibit- 
ing sources due to non-invariant transformations can be 
designed using electric charges and currents, magnetic 
charges and currents or both, which clearly character- 
izes physically different situations with the same media 
properties x- 

Does there exist the possibility of a geometric interpre- 
tation of Eq. ([75]) ? If this shall be possible space must be 
transformed differently for different components of F^'^ 
and -ffpj/. In fact, the most general linear transforma- 
tion can be interpreted as a separate transformation of 
spacetime for each component of the two tensors. Let us 
provide a simplified example, where independent spatial 
transformations are applied to E, B, V and Ti. Labora- 
tory space is denoted by a;\ the electromagnetic spaces 
by a;^, a;^, x^j-, and x^^j, resp. Under time- independent 
spatial transformations all four fields transform as (co- 
)vectors and thus Eqs. (|A.9p - (jA.12p suggest the inter- 
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pretations 



Ei = se „ .- E^ 



dx^E ^ ' 



B* = 



V 



7d dx 



D 



y^ dx^ 



IB dx^g 

-^7 dx^ 
~ dx-^ 



B^ 



dx: 



H 



yielding for permittivity and permeability 



e*^ = SE- 



/TD dx'-o j,i dx'j^ 



y^ dx^ dx^ 






/7 dx^ 



dx' 



(81) 
(82) 

(83) 
(84) 



As is seen from ([1]) and ([2]) Gauss' law for B^ and 2?* 
remain unchanged, while Faraday's and Ampere's laws 
are changed according to 



IE dx 



IB dxi 



^Wo&+e''''d,Ek = 0, 



,^^'^O,n,-^^VoV^=0 



/iD dxi 



D 



making the electric and magnetic currents 



VoV' 




(85) 



(86) 



(87) 



necessary. For completeness it should be mentioned that 
the transformations ([5T1) and ([5^ allow a straightforward 
interpretation since each of the four Maxwell's equations 
still can be transformed as a whole. Taking even more 
general transformations, e.g. transforming each compo- 
nent of the electric and magnetic fields separately, does 
no longer allow this manipulation in a simple way and 
thus will make the derivation of the necessary media pa- 
rameters more complicated. 



light rays travel along trajectories as if the medium was 
a transformed, empty space. Still, the transformation 
medium in general is quite different from transformed 
empty space, if the conservation laws from the stress- 
energy-momentum tensor are considered. This aspect is 
even more important within the extension proposed here, 
as there exist two different transformed (electromagnetic) 
spaces and light rays don't follow the geodesies of any of 
them. We have shown that one can make a virtue out 
of necessity: the geometric approach does not just pro- 
vide a tool to design the path of light in a medium, but 
equally well it may be used to design the behavior of 
(parts of) the stress-energy-momentum tensor, e.g. the 
direction of the Poynting vector, and/or the behavior of 
the wave vector. Here the proposed generalization of- 
fers many more possibilities compared to the known dif- 
feomorphism transforming media. In particular we have 
derived the geometric relations that describe the trans- 
formation of the Poynting vector, of the direction of the 
wave vector as well as the dispersion relation. 

Finally we have commented on a different route to 
generalize the notion of transformation media [lOJ, llll |. 
These field-transforming media are not based on invari- 
ant transformations of the equations of motion and con- 
sequently source free solutions of the original configura- 
tion are not mapped onto source free solutions of the 
new medium. We have shown that also this approach 
may be covered by a generalized concept of coordinate 
transformations. Still, there remains a fundamental dif- 
ference between the approach of Refs. [10, 11] and the 
one discussed here: While in the former case the trans- 
formations are ultra-local (the transformed fields at the 
point x^ are defined in terms of the original fields at this 
point), in the latter they are essentially non-local, as the 
transformed fields at x^ are related to the original fields 
at some x^ ^ x^ . The preferable approach depends on 
the specific problem at hand, also a combination of the 
two is conceivable. 



VII. CONCLUSIONS 



In this paper we have introduced a generalization of 
the concept of diffeomorphism transforming media, the 
basis of transformation optics 1, 2, i3i|. As basic idea we 
have found that spacetime can be transformed differently 
for the field strength tensor (containing E and i?) and the 
excitation tensor (encompassing T) and 7i). This exten- 
sion allows design of non-reciprocal media, in particular 
permittivity and permeability need not longer be sym- 
metric. Furthermore, this approach permits a geometric 
interpretation of indefinite media [g, |7| . 

Diffeomorphism transforming media are motivated by 
the wish to produce a medium that looks like a trans- 
formed but empty space. The basis of this interpreta- 
tion is Fermat's principle applied to these media 'j||: in- 
deed it is found that in a transformation medium the 



Acknowledgments 

The author wishes to thank J. Llorens Montolio for 
helpful discussions. This work profited a lot from fruitful 
discussion with C. Simovski, S.A. Tretyakov, I.S. Neve- 
dov, P. Alitalo, M. Qiu and M. Yan during a cooperation 
of the Advanced Concepts Team of the European Space 
Agency with the Helsinki University of Technology and 
the Royal Institute of Technology (KTH). The coopera- 
tion was funded under the Ariadna program of ESA. 



APPENDIX: COVARIANT FORMULATION 

In this Appendix we present our notations and conven- 
tions regarding the covariant formulation of Maxwell's 
equations on a possibly curved manifold. For a detailed 
introduction to the topic we refer to the relevant litera- 
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ture, e.g. [T3,[i3|- Throughout the whole paper natural 
units with cq = ^q = c = 1 are used. 

Greek indices ^^v, p, . . . are spacetime indices and run 
from to 3, Latin indices i, j, fc, . . . space indices with 
values from 1 to 3. For the metric we use the "mostly 
plus" convention, so the standard flat metric is 5^1/ = 
diag(— 1, 1,1,1). If we interpret x^ =t with (laboratory) 
time the space metric can be obtained as \\^ 



tory space 

dx^ dx^ . 

^ ax^ ax^ ax'' 



1^' = 9^' 



li] = 9lk 



goo 



I'^ljk = SI 



V 



g / / dx^ dx'- dx° dx^ \ 



(A.l) 

This implies as relation between the determinant of the 
spacetime metric, g, and the one of the space metric, 7, 

-9 = -goal (A.2) 

The four dimensional Levi-Civita tensor is defined as 

1 



— 5 \ \ dx^ dxi dxi dx^ J 



X)j 



dx° dx^ 
dxi dx'' 



T^'^'ni I , 



(A.9) 
(A.IO) 
(A.ll) 



^ == V-900 dxJ 
Hi = sa -^^Hj 

y/-90Q ox"- 



g dx^ dx'' I 



/^ dx^ dx' 



We characterize the general linear, lossless media usu- 
ally by means of the Tellegen relations 



e^.pa = V^[fiiyp<j] , et'-'P^ = — ^fiiypa] , (A.3) V ^ e'^E, + k'^j , B' = p'^Hj + C' E, . (A.13) 



with [0123] ~ 1. Therefore the reduction of the four 
dimensional to the three dimensional tensor reads 



In terms of field strength and excitation tensor the media 
relations become the Boys-Post relation 



^Qijk — \/ — 900^ijk , 



Oijk 



1 



H"" 



1 



A ^ pry y 



(A.14) 



ijk 



V-900 



(A.4) 



where x niust be invertible with inverse 



An additional complication arises in the definition of 6^^. 
and eijk, since the orientation of the spacetime manifold 
may change without changing the orientation of space 
(e.g. by a mapping t = —t.) Therefore the corresponding 
relations should be written as 



Xf^.XrX^^'"' = {S'pS: ~ S;SP) . (A. 15) 



By virtue of Eqs. (|A.6p and (|A.7p Eq. (|A.14p may be 
written as 



soijk = (^V—goo^ijk , eoijfe — cTy—goo^ijk , (A. 5) 

where tr = +1 if space and spacetime have the same 
orientation and a — ^1 otherwise. 

The field strength tensor F^j, encompasses the electric 
field and the magnetic induction, the excitation tensor 
H^'^ the displacement vector and the magnetic field with 
the identification: 



Hi 
-pi 






Ej 
Bi 



with 



21'^ = -^/^9ooX 
1 



OiOj 



*Sij — -\/ —goO^ikl^jmnX 

U' = -\^^^o^^klX^'''' 



Ei = Foi , 



B' = --e'^'^F. 
2 



jk 



(A.6) 



s^ = -^^-goo^jkix 



Oikl 



(A.16) 

(A.17) 
(A.18) 

(A.19) 

(A.20) 



■Hi 



\/-9oo 



ei.kH"' . (A.7) 



The Tellegen and Boys-Post formulations are related by 

e'^ = (a - d^-^g:)*^ te^ = (ds^^)'^ (A.21) 



Finally, electric charge and current are combined into a 
four-current J^' = {p/^/-goo,f/V-9oo)- F^,^ and H^"' 
are tensors, thus under the transformations of Section Hill 
they behave as 



p'^ = (» 



-ivj' 






{%''€)'' 



(A.22) 



F„, 



dx^ dx" 

d^ ""ai^ 



JJP^ 



dx^jjpa^ 

dxP dx° 



(A., 



This implies for the transformed space vectors in labora- 



Finally, we mention that the rank 4 tensor x''"'^'"^ may be 
decomposed as [l^ ]1^ 

p,upa _ (1) ui/pcr j^ ^pu\[p g a] _ ^pa\[p g v] _|_ Q,gM'^PO- 

(A.23) 

where the principal part '^^'x'^'^'"^ has no part completely 
anti-symmetric in its indices and is symmetric under the 
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exchange (/i, i^) <-^- {p,<7). The principal part has been 
discussed extensively in [l2|. 5'^'' was introduced in 
Refs. [il, [I^ as skewon part (related to chiral proper- 



ties of the material 1161 ,1211 ) , while a represents the well- 
known axion coupling [22| . 
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[^ differs from the one used here. According to our con- 
ventions F^" = Z/'^" in vacuum, while there ^J—gF^" — 

Sources are external parameters and thus should be set to 
zero for a symmetry transformation. Even together with 
sources the symmetry can be restored, if an appropriate 
transformation rule of the sources is defined. 
Notice that under the continuous transformation the ac- 
tion behaves as 5 ^ (cos^ a — sin^ Oi)S and thus for 
a = 7r/4 transforms to zero. Therefore at the level of 
the action only the discrete duality transformation can 
be considered. 

"Time" t here is just a variable to parametrize the world 
line, this choice does not make any statements about the 
speed of light. 
Alternatively, one could start from a relation of the form 

//''" = A^^p^H"" + B^^P^Fp^ 

but the form H76p appears more transparent to us. 



